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Abstract 

We argue that supersymmetric transformation could be applied across the heterojunc- 
tion formed by joining of two mixed semiconductors. A general framework is described 
by specifying the structure of ladder operators at the junction for making quantitative es- 
timation of physical quantities. For a particular heterojunction device, we show that an 
exponential grading inside a nanoscale doped layer is amenable to exact analytical treat- 
ment for a class of potentials distorted by the junctions through the solutions of transformed 
Morse- Type potentials. 

PACSnumber(s): 03.65.-w,03.65.Ca,0.65.Ge,73.40.Kp,73.40.GK,73.40.Ty,73.43.Jn,02.30.Hq,02.30.Gp 

To have a proper grasp of designing any semiconductor device that relies on a heterojunction, 
a complete understanding of the transport phenomena at interface is important [1] . Positive and 
negative carriers moving across a junction [5J [3] face an abrupt jump in potential due to the 
periodicity defects in crystal lattice at interface and varying electrical and optical properties 
of different materials bounding the heterojunction. Different values of the position-dependent 
conduction-band energy on the two sides separated by a junction results in a <5-likc singularity 
at the interface j4j producing either a very deep, narrow well or a very high but a thin barrier 
that depends on the position of the junction. In these situations it is worth considering the 
transformed potential which in a supersymmetric context shares almost the same bound spectra 
but are free from such a singularity. A general version of supersymmetric transformation (see the 
monograph [5] for an up-to-date survey) allows to construct a hierarchy of potential chain which, 
in physical purpose, is useful to remove or embed arbitrary number of bound states. Since the 
position of a bound state is related to the pole-structure of the S matrix, its artificial designing 
can be utilized to improve the mobility of the carriers and thus having a desired output at the 
acceptor diode. A study of the spectral properties of the induced intertwined potentials across 
the heterojunction is thus of immense physical importance. 

In actual microscopic description of the motion of electrons (or holes with e — > — e) across 
a semiconductor junction, it is known [6] that the supersymmetric structure of the underlying 
Dirac Hamiltonian describes unpaired 'spin-up' and 'spin-down' states in the conduction and 
valence band near T or L-point in the Brillouin zone that are localized at the junction. Note that 
the existing literature [7-11] does not apply to a practical semiconductor device where two graded 
mixed semiconductors are artificially glued (say Gai_ x Al x As/GaAs heterojunction or Si-Si02 
layer) since the role of ladder operators across the junction is unclear. Consequently the usual 
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techniques of constructing the transformed potential fail in the effective-mass approximation 
which replaces actual wave function by envelope wave function. In this Report we show that the 
information about spectral properties of a given effective mass Hamiltonian having a (5-singularity 
at the hetcrojunction can be extracted from the supersymmetrically transformed Hamiltonian 
which is free from such a singularity. We apply our method for a particular Gai-xAl^As/GaAs 
semiconductor device to compute analytically local electron density based on inter subband 
energy levels and envelope wave functions obtained from a physically justifiable choice of grading 
function. 

We start by considering two heterojunctions formed by an artificial joining of three different 
material layers (doped or pure). In the effective mass approximation the normalized envelope 
wave function of a single electron in the n-th subband moving along the normal to the surface 
of the material layer (xy-plane) obeys the Schrodinger equation |12| for the given Hamiltonian 
Hq(z) in SI unit 



2m* (z) dz 



+ V {z))i) n {z) = E n tp n {z) (1) 



The effective mass m* has a position-dependent variation mi(z) between two junctions at z = 
ai,ci2 connecting continuously two known constant values m,\ and vn^ outside the junctions 
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m*(z) = Gjrrij ; TOa(ffli) = mi , i = 1,2. (2) 
In above equations GjS are dimcnsionlcss Hcaviside theta functions defined by the relations 

Gl = Q (^£) 5 G3 = e(i^a ); G2 = i-(g 1 + g 3 ), (3) 

where a = 4:ire h 2 /moe 2 = 0.0529 nm is Bohr radius, e denoting permittivity of free space and 
mo,e are free electron mass and charge. A general class of the potential Vq(z) containing the 
(^-singularity at the junctions is 



/ _ dm 2 i 

m 2 = — ,e, = (-l) 



(4) 



The actual form of the component potentials Vq(z) will depend on chemical composition of 
materials in different regions and the applied electric and magnetic field. 

We assume that inside each region (except for the junctions) the component Hamiltonians 
= -d/dz[(h 2 /2m 3 )d/dz] + V£ j) possess formal solutions ipjp (z) for an arbitrary energy E. 
We shall prove that there exists transformed Hamiltonian Ho(z) that shares the same spectra 
as Hq but is free from ^-singularity. In other words using supersymmetric factorization we are 
going to construct a new potential of the form 

Vo(z) = Y,GjV ij \z), (5) 

~U) 

satisfying the above properties where Vo {z) are non-singular. Note that the component Hamil- 
tonians can always be projected in a factorizcd form 
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where the superpotential W is a continuous smooth function determined by the shape of the 
component potentials. Our aim is to factorize the given Hamiltonian Hq in terms of the ladder 
operators which will act globally through the junction, i.e. Hq = AA^ — £ (factorization energy 
£ is chosen to make spectra positive definite). 

Since all junctions are quantum mechanically of equal preference it is reasonable to choose 
the form of the ladder operators at each junction as the average of the two adjacent component 
operators. This leads to the following representation of A: 

A = YG 3 A ] = ^=^ + W{z) (7) 

It is straightforward to verify that the presence of theta functions in the ladder operators simu- 
lates (5-singularity in the given potential Vq(z). It remains to show that the transformed Hamil- 
tonian given by the factorized form 

H {z)=A*A-E (8) 

will not contain such a singularity. Clearly the wave function ip of the given Hamiltonian 
must be continuous in the entire region for otherwise Schrodinger equation ([1} has to hold a 
stronger singularity than 5. Noting that the intertwining relation 

H A = AHo (9) 

relates ip and ip as ip oc Aip, we conclude that the action of A on ip must be continuous across 

junctions. Suppose that the the wave functions ip of Hq may be constructed as a superposition 

7(i) 

of component solutions ip^ where the average of two adjacent plain waves contributes equally 
at the junction. Using this construction and the representation (0 of the operator A we then 
arrive at the following matching conditions on ip 

$'(0j + 0) = ${aj - 0) = i)'{aj). (10) 

where we have used the usual continuity condition for ip. The smoothness of ip across junction 
indicates that transformed potential (|5|) is free from any singularity. This can be directly verified 
from the explicit expression of the transformed potential 

Vo(z) = w*-(^Ly-z. (ii) 

If one applies similar reasoning for ip oc A'lp, a slope-discontinuity of ip will appear at the 
junctions which is the signature of ^-singularity. This striking difference between two transformed 
Hamiltonians arises due to the non-commutativity of momentum and mass operators which we 
can exploit via supcrsymmctric transformation. 

Equations (|7|.(|10p and (fTTj) constitute the central result of this Report establishing a way 
to apply supcrsymmctric transformation across heterojunction. The set of conditions in (|10j) 
yield a transcendental energy equation whose roots correspond to the allowed electronic energy 
levels for the given Hamiltonian except for the factorization energy £. Finally we note that for 
practical application of above result the particular composition of material layers will guide the 
choice of the superpotential simulating the potential structure inside and outside the junction. 

Let us consider a device made up of Gai-^Al^As/GaAs with the interface located at z = 0. 
The n-type Gai-^ALj layer is grown in two Regions I(— Zh < z < —at) and II(— at < z < a t ) 
and adjacent GaAs layer is grown in Region III(at < z < z t ) so that two heterojunctions are 
symmetrically placed about the origin. The potential Vq(z) is generally formed from several 
contributions 

V {z) = -ecj>{z) + V h (z) + Vi m (z) + V xc (z) (12) 
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The electrostatic potential <f>{z) satisfies Poisson's equation (SI unit) 



d_ 

dz 



e K(z) — 
dz 



en(z), (13) 



where n(z) = J2 n ^n^Pn( z ) ~ Nd(z) + Nac(z) is the electron concentration in which 
Ne>(z),Nac{z) denote the position-dependent donor and acceptor concentrations and N n (z) 
represents the number of electrons per unit of area (in units of cm _2 )in subband n: N n (z) = 
{m* {z)ksT / ttIt 2 ) ln[l + exp{(£i? — E n )/ksT}, Ep denoting Fermi energy. Since electrostatic 
potential has to be extrapolated between ([1]) and (|T3"| in an iterated numeric procedure, we shall 
not consider it here. Vi m (z) and V xc (z) are respectively image potential [arising due to the dielec- 
tric constant step k{z)] and exchange-correlation potential due to electron-electron interaction 
in the channel. It is well-known that image potential cannot be incorporated in equation ([1} as 
it involves singularity at the origin. This effect is usually realized into V xc (z) by adding an image 
term dependent on mutual coordinates of two interacting electrons in the population. However 
for Kb ~ k c [the subscripts b and c indicate the barrier (Gai-^Al^As) and the channel (GaAs) 
sides of hetero junctions], image term could be ignored and exchange-correlation potential may 
be parameterized [T3] as (in Hartree atomic unit of energy: Eh = ft 2 /moap ~ 27.2 eV) 

, , 1.,,, .2£m*(z)/m 
V m (z) = -^K 2 cr {x) E h , (14) 

where the dimensionless constant K cr (x) depends on the doping fraction x = x(z) as 




0.7734xln(l + x-i) 

21^(4/9^/3 ( 15 ) 

For the analytical purpose, a suitable continuous functional form for n(z) is useful which con- 
nects two known values Kb and k c of dielectric constant. A simple choice k(z) = \J '2£m* (z) / mo 
reduces V xc {z) to an x-dependent constant simulating doping-innuenced shift in the energy scale 
where £ is a parameter adjusting with Kb and k c . Usually in numerical process a linear interpo- 
lation is taken between the values in the barrier and channel side. We propose an exponential 
grading in the intermediate Region II to connect values in I and III. This is reasonable since in a 
sufficiently narrow thickness of layer II (2at in unit of nm) it basically implies linear interpolation 
and add corrections with the increase of thickness. Hence following connection between known 
effective mass values mb = rn2(—a t ) = brriQ and m c = m2(a t ) = cmo on Region I and III will be 
adopted 

m*(z) = G 1 m b + G 2 m2 + G 3 m c , m 2 = m p 2 e 2Pz/a ° (16) 

where the dimensionless parameter /3 is related with mass-steps b, c on both sides and the grading 
constants are given by ([3]) with the replacement of ai,2 by ^fat- The choice of the supcrpotcntial 
will be guided by the physical fact that Aluminium doping creates a barrier |14j for electrons 
to diffuse to Region I and so they are mobilized to acceptor site containing a GaAs quantum 
well. Thus the potential Vq should have a barrier on left, a well at intermediate region and 
approaching a constant value on the right to make the electrons pass through. This situation 
can be realized by the following supcrpotcntial 

W(z) = ^=(A - Se-^/ ao ). (17) 

a V2m 



Both parameters A and 23 are dimensionless. The latter parameter is free controlling the 
depth of well and the height of the barrier and the former is related with exchange-correlation 
term fj 14[) through the definition A = K 2 — K 2 r , K > K cr . At critical value A'= K cr , the 
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Figure 1: Original and transformed potentials with three levels Eq — 43.58, Ei = 46.75 and E2 = 51.67 
eV (E F = 35.9 meV at 300 K) for the choice A = 2, "B = 1.1, /3 = 0.21 and a t = 4.73 which corresponds 
to 0.5nm intermediate layer thickness. Electron concentration is calculated using the data indicated in 
figure based on three wave functions obtained. Electrostatic term is not included. 



doping-influence will be prominent, since it has very small values for x = 0.3 and 0.4. Our 
strategy is to solve the transformed potential ([3]) according to the matching condition (|TU1) . 

The add-on advantage of the model-choice (|16[) and (|17[) is that the component potentials 

~(i) . 

Vo in three regions are of Morse type and so exact analytical solutions are available. The final 
expressions of unnormalizcd solutions for the transformed Hamiltonian are 

We \ z ) = exp[-^/b{x(z) + ^/m Q Xz}]U(g b ,h b ,2^/bx) 

^\z) = exp[^-y(z)][M(g 2 ,h 2 ,2y(z) + gU(g 2 ,h 2 ,2y(z))} I (18) 
4>e\z) = exp[- v / c{x(2) + y/m Xz}]M(g c ,h c ,2y/cx) J 

where X {z) = T,e- 0z / a ° / ft, y(z) = a X y /m e l3 ^ ao , X(E) = y/[A 2 /m al - 2{E + £)/h 2 ] and 
U(g,h,£), M(g,h,£) are two linearly independent solutions [TS] of confluent hypergeometric 
equations 

of -—o m 

with the arguments <7j = ^/j(ao^moX—A)//3,hj = 1 +2aaXJm~j /j3;j — b,candg 2 = (223 — 1)(A— 
A)/2X, h 2 = 223 — 1. The coefficient g in the second equation of (|T5|) is evaluated at solutions 
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of energy equation. The normalized envelope wave function of original Hamiltonian (HJ for n-th 
subband E n (> 0) is 



3=1 



(20) 



w 

a 
o 

■H 

+J 


a 
p 

h 

CD 
> 

s 

CD 

N 
-H 

H 

td 





4125 



1375 



-1375 



-4125 



-6875 



<-Region I-> -Regi 

/ 

/ 


on II -|| 

^2 

\/ A 


<- Region III -> 

/ X 3 










A I 








\ // Supersymmetric Transformation 






i// n =Ai// n+1 , n=0,l,2 


- i i 







-0.79 



-0.25 0.25 0.79 1.32 

Distance z (nm) 



1.85 



Figure 2: Envelope wave functions for the original and transformed potentials. Wave function for given 
potential is obtained by applying ladder operator A which suppresses one node. Electrostatic term is 
not included. 



The energy equation is solved, using Mathematica, for representative values: A = 2, 23 = 1.1 
and p = 0.21 with the intermediate layer thickness 0.5nm. Fig [T] shows the given potential 
(solid line) with <5-singularity and transformed potential (dashed line) with J-singularity removed 
having four levels with zero energy ground state. Thus we choose £ = and the three levels 
for given potential are shown. The levels are very close (E\q = 3.17, E21 = 4.92 eV) and lie at 
the bottom of the well above the Fermi energy Ep = 35.9 meV at 300K. Electron concentration 
(grey curve) is shown at 300K by choosing smooth continuous functional form around small 
intervals containing two junctions. Fig [2] shows the effect of supersymmetric transformations on 
wave functions Vn+i> n = 0, 1, 2 (Grey curves) which suppress one node to recover original wave 
functions ip n , as is expected. Higher levels for transformed potential violate oscillation theorem 
as ^-singularities resist creation of node near junctions. Interestingly original wave functions 
perfectly follow oscillation theorem leaving a sharp non-smoothness at the junction near channel 
side. Hence transformed Hamiltonian, which is free from ^-singularity, digests its effects and 
makes the original solutions physically acceptable. 
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To conclude, we have provided a method to bypass the essential (^-singularity in the effective- 
mass potential at the heterojunction using supersymmetric transformation. Note that the 
method is initiated from the observation that although original potential is invariably affected by 
(^-singularity at the heterojunction, its supersymmetric partner escapes such a singularity. We 
emphasize that this observation was unnoticed in the literature. Further, in spite of extensive 
works [7-11] regarding application of supersymmetry in position-dependent mass problems, no 
prescription existed to tackle the presence of heterojunction [5] in semiconductor. This gap is 
filled in the present communication. Finally, for a particular Gai^Al^As/GaAs device we have 
provided an exact analytical treatment based on a model grading function connecting values 
of band parameters across the junction. The choice of exponential grading function is subtle 
for exact analytical solution as it leads to a mapping in a constant-mass scenario to a known 
Morse potential in Region I and III which are joined by a distorted Morse- Type potential in the 
intermediate Region II. 

We thank Dr. S. Ghosh for his interest in this work. AS is supported by DST, Govt, of India. 
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